# -*- coding:utf-8 -*-

def calc_wpd(ws, rhoair=1.125):
    return 0.5 * rhoair * ws * ws * ws

def calc_rhoair(p, t):
    """
    :param p: Pa
    :param t: K
    :return: kg/m^3
    """
    rhoair = p / (287.058 * t)
    return rhoair


